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This paper makes two points. First, we show that the line-of-sight solution to cosmic microwave 
anisotropies in Fourier space, even though formally defined for arbitrarily large wavelengths, leads 
to position-space solutions which only depend on the sources of anisotropies inside the past light- 
cone of the observer. This foretold manifestation of causality in position (real) space happens order 
f^ by order in a series expansion in powers of the visibility 7 — e~'^, where fi is the optical depth 

y—i to Thompson scattering. We show that the contributions of order 7^ to the CMB anisotropies 

^^ are regulated by spacetime window functions which have support only inside the past light-cone of 

Cn the point of observation. Second, we show that the Fourier-Bessel expansion of the physical fields 

j^.,^ (including the temperature and polarization momenta) is an alternative to the usual Fourier basis 

jrt as a framework to compute the anisotropies. The viability of the Fourier-Bessel series for treating 

*^ the CMB is a consequence of the fact that the visibility function becomes exponentially small at 

<H redshifts z ^ 10^, effectively cutting off the past light-cone and introducing a finite radius inside 

.^^ which initial conditions can affect physical observables measured at our position x = and time io- 

Hence, for each multipole i there is a discrete tower of momenta ki i (not a continuum) which can 
affect physical observables, with the smallest momenta being kn ~ £. The Fourier-Bessel modes 
take into account precisely the information from the sources of anisotropies that propagates from the 
initial value surface to the point of observation - no more, no less. We also show that the physical 
observables (the temperature and polarization maps), and hence the angular power spectra, are 
f^ unaffected by that choice of basis. This implies that the Fourier-Bessel expansion is the optimal 

^L( scheme with which one can compute CMB anisotropies. 

6 

I 1 I. INTRODUCTION 

►;>. The cosmic microwave background (CMB) is the earliest, cleanest observation that reveals what the Universe looked 

^f) like at the very beginning. A remarkable string of observations of the CMB temperature fluctuations over the last 
^O 20 years, most notably by COBE-DMR [1] and WMAP |2}jJ,, has shown that the typical initial conditions of the 

'^ Universe when it was under 400.000 years old can be characterized by extreme homogeneity and isotropy, only slightly 
^— ^ perturbed by small, 0(10^^) fluctuations with a nearly scale-invariant spectrum. More recently, the small degree of 
l/~i polarization that is imprinted on the CMB radiation by anisotropic Thompson scattering has also started to become 
(^^ detectable [5ll7], and may hold the key to unravel the mistery of the birth of our Universe - see, e.g., JSHTT]. 
^^ The remarkable success of the CMB as probably the most powerful tool in observational Cosmology can also be 

, . traced to the simplicity of its underlying mechanisms: Thompson scattering and linear perturbation theory. The 
^ basic theory, which is an application of the relativistic radiative transfer equations |12| , was initially developed in 
connection with the CMB by Peebles and Yu [13] , and the flrst to write down the full coUisional Boltzmann equations 
for temperature and polarization were Bond and Efstathiou [111 I15j. 
j^ However, even if the main mechanisms driving acoustic oscillations in the baryon-photon fluid were basically under- 

stood early on, crucial features such as neutrinos, gravitational waves, spatial curvature and the effects of lensing on 
polarization remained puzzling. It was not until the 1990's that the theory reached full maturity, when the complete 
separation of radial and angular modes allowed by the use of spin angular momentum eigenfunctions for polarization 
cleared the way for our current understanding of CMB physics |16fBl] . 

One particularly signiflcant step forward was achieved with the line-of-sight solution to the coUisional Boltzmann 
equations |25j. The idea is that photons travel along null geodesies, hence the comoving distance Ax between two 
successive collisions is equal to the conformal time interval Ary between those collisions (as usual, we assume that 
the Born approximation is valid.) This means that, given a line-of-sight h, we know that a photon detected at time 
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r] travelling along that line-of-sight was at the position x' = nirj — -q') at time rj' , if no collisions occurred between 
those times. The more general case of an ensemble of photons can be easily accomodated in this picture, since the 
probability that a photon scatters with free electrons is given in terms of the optical depth for Thompson scattering, 
which to a very good approximation is a smooth function of time, ^{ri). The final state of the ensemble is then 
obtained through integration over time of the sources of temperature and polarization anisotropics, appropriately 
weighted by the optical depth at each time. 

The power of the line-of-sight integral solution is that it separates, as much as it is possible, the (free) propagation 
of photons from the ultimate sources of anisotropics (matter and metric perturbations) - so, it is similar in spirit to 
a Green's function for the temperature and polarization of an ensemble of photons. However, there is one feature 
of the generation of anisotropics which makes it impossible to completely separate the sources and the anisotropics: 
anisotropic Thompson scattering is itself a source of both temperature and polarization, so the process is, in some 
sense, non-local. 

When photons scatter at a given place and time, the fluctuations in temperature and polarization that are generated 
as a consequence of those scatterings depend also on the quadrupole of the temperature and the polarization of the 
photons that were incident at that place and time. Therefore, the fact that those incident photons typically propagated 
to the location of the scattering from far away implies that the process is non-local - hence the line-of-sight solution 
is actually a set of integral equations, at least for the lowest multipoles {i < 2). The higher multipoles {i > 3), 
however, can be completely determined from the lowest multipoles, which makes the line-of-sight solution a vastly 
superior method compared to the usual hierarchy of Boltzmann equations. Of course, for the low multipoles the 
integral equations are impractical, and the preferred method to compute them is to revert back to the hierarchy of 
Boltzmann equations, which is then truncated at a relatively low multipole that is sufficient to accurately compute 
the multipoles £ < 2. 

Nevertheless, despite the fact that the generation of anisotropics is a non-local mechanism, it is still completely 
causal: photons propagate along light cones between scatterings (which are basically instantaneous within the cos- 
mological timescales.) The main drive behind this work is to clarify how causality and non-locality are manifested 
in the generation of the CMB. We are interested, in particular, in describing the generation of the CMB in position 
(real) space, since Fourier space can sometimes obfuscate the causal nature of the physical mechanisms. 

We will show that, in position space, the non-local anisotropics can be resolved into explicitly causal pieces by 
iterating the integral line-of-sight equations. The smallness parameter in this perturbative expansion (which is rem- 
iniscent of the Dyson series of Quantum Electrodynamics |i26j ) is the total visibility 7 = e~'^, where ^ is the optical 
depth to Thompson scattering. The first term in that series corresponds to the last scattering of the photons before 
they were detected; the second term corresponds to the last two scatterings before detection; and so on. 

In this series over the number of scatterings, causality is manifested at each order in terms of radial integrals corre- 
sponding to spacetime window functions which are only non-zero inside the past light-cone (PLC) of the observation 
point. The spacetime dependence of the 7^-order term is given in terms of integrals of iV-l-l spherical Bessel functions, 
which we demonstrate to have support only inside the PLC. The first term of the series, 0(7^), corresponds simply 
to the Sachs- Wolfe (SW), Doppler and Integrated Sachs- Wolfe (ISW) contributions to the temperature fluctuations 
from the last scattering [27J, and evidently it is only nonzero at the surface of the light cone. The second-order terms 
for temperature and polarization carry the memory of the last two scatterings, and can be non-zero anywhere inside 
the volume of the PLC (not only on its surface). Fig. fl] illustrates the structure of the light cones for one, two and 
three scatterings. Hence, we have shown that the line-of-sight integrals are actually retarded Green's functions for the 
temperature and polarization, and we found analytical expressions for them in position space. 

The second result of this paper follows from this causal structure. The expansion in successive scatterings (or, 
equivalently, in orders of the visibility function) makes it clear that what matters for the CMB observables are not 
simply the source fields and the anisotropics as functions of position and time, /(x, 77), but the fields multiplied by 
powers of the visibility, jirj)^ f{x, 77). Since the visibility vanishes for very early times (say, Zi >> 10'^), for all practical 
purposes that have to do with the CMB, our PLC is cut-off at rji — 7]{zi), and all the sources are effectively zero 
outside the radius corresponding to that time (which in a typical flat ACDM model is approximately r^ sa 5Hq .) 

The fact that the physical fields are exponentially suppressed at the boundaries of the PLC means that the best 

basis for expanding the fields, as well as the anisotropics, is not the Fourier basis, which is most convenient for plane 

waves in M^, but the Fourier-Bessel basis [551, which expands functions f{x) inside a 3D sphere of finite radius R into 

spherical harmonics Ygm{x) and spherical Bessel functions ji{kiix), with k^R being the i*'' root of j^. Both sets of 

-a „ 

basis functions are eigenvectors of the Laplacian operator in fiat space, with eigenvalues —k or ~k~^. but they differ 

most significantly in that the Fourier-Bessel series establishes a discrete tower of momenta for each multipole I, so that 

the smallest momentum mode is ku ~ t/Ti. Therefore, there is a clearly defined minimal mode that can contribute 

to CMB observables, and the modes above it are all discretized. The drawback is that, if the underlying fields are 

Gaussian in nature, the Fourier-Bessel modes, as opposed to the Fourier modes, are not statistically independent, 

so their covariance matrix is not diagonal. However, dynamics is one thing, statistics is another, and we will show 
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Figure 1: Past light-cones for photons experiencing one, two and three scatterings since the time of decouphng from matter 
(denoted as the base of the cone.) The terms corresponding to these interactions are respectively of order 7^, 7^ and 7^, where 
7 — e~^ is the visibility. 

how these issues can be separated so that we can easily recover angular power spectra which have precisely the same 
statistical properties of the power spectra computed in Fourier space. 

The first part of this work takes an approach which is similar to that used to make constrained simulations of 
CMB temperature and polarization maps by Liguori, Matarrese and Moscardini |29| and by Komatsu, Spergel and 
Wandelt [3^ - see also |31fl55] . However, while those simulations make use of several approximation schemes in order 
to produce numerically viable codes, our expressions are exact and analytical. In particular, our results show that 
the transfer functions of [29, "31? have an invariant (or geometrical) piece which is described by our spacetime window 
functions, so that the transfer functions can be obtained by integration of these window functions over time with 
some visibility function. Another work close in spirit to ours was done by Bashinsky and Bertshinger |34l 135] , who 
calculated the Green's function for the evolution of linear cosmological perturbations in position space but did not 
compute the Green's functions for the anisotropics. 

A quick note on our conventions: the Fourier transform of a function f{x) is f{k) = {2tt)^^'^ J (fix e"'''''^ f{x): 
the perturbed Friedman-Robertson- Walker metric is ds^ = a'^{i])[—{l + 2^)dri^ + (1 — 2'i)dx^]; and the relationship 
between the rotation matrices and the spin spherical harmonics is D^j _g{a,(3,j — > 0) = •y/47r/(2£ -I- 1) sY^^{/3,a), 
so that sY^m ~ (—1)'"^" ^sYe.-m- Finally, everything we will say in this paper concerns scalar (spin-0) density 
perturbations - we will consider lensing and gravitational waves in future work. 

II. CMB ANISOTROPIES IN A NUTSHELL 

The main channel of electron-photon interactions during recombination is elastic Thompson scattering. The likeli- 
hood of a photon interacting with free electrons between times 771 and 772 can be determined in terms of the optical 
depth for Thompson scattering: 

rV2 
m('7i,'72)=/ dr]a{ri)aTne{ri)Xe{r]) , (1) 

Jrii 

where a{rj) is the Friedman-Robertson- Walker scale factor, ar — 6.65 x 10^^^ cm^ is the Thompson cross section, Ue 
is the total number density of electrons, and X,, is the ionized fraction. The probability per unit time that a photon 
observed at 772 had last interacted at time 771 is called the visibility function: 

5(^,^2) = /-e-^^''^^''^^ (2) 



where g{r]i,r]2) is a positive-definite function, normalized to unity if the limits are taken such that /i — >■ and /i — ^ oo 
at late and early times, respectively. It is customary to define the visibility function today simply by g{r]) = g{r],rio). 
In this work we will also define the total visibility (or unscattered fraction) as: 

7(^1,^2) ^e-^^"-"^). (3) 

Notice that the total visibility is the probability that a photon will not scatter after rji before it is detected at r/2, and 
in fact we defined it such that it is related to the visibility function by g(j]i,r]2) = d/dr^i [7(771, 772)]. 

At very early times, the optical depth is extremely large, so the Universe is effectively opaque and the visibility 
vanishes exponentially. Photons and electrons are interacting so often that the photons can spend enough time in 
some small region so that inelastic processes lead to thermalization with the electrons and, by extension, with the 
baryons. Hence, at very early times (much before recombination) the photon distribution @(x,ri;l) = AT/T was 
essentially in equilibrium with baryonic matter, and only the monopole 0Q(x,r]) = J (Pl/{4Tr)Q{x,ri;l) = 5j/A was 
significant (here 6^ is the density contrast of photons.) 

However, as soon as recombination starts the optical depth plummets, scatterings between photons and electrons 
become sparser, and as a result the radiation incident on any given scattering source can be increasingly anisotropic, 
since photons arriving to a scattering source from distant over- or underdense regions have different equilibrium 
temperatures. Thomson scattering with anisotropic radiation then generates polarization (and vice- versa), and the 
process becomes quite intricated. 

Let us assume for a moment that the photons decoupled from matter instantly, at some time 77^, and propagated 
freely from that time down to our detectors at x = 0, 770 . Then direct integration of the geodesic equation (with 
the help of the Born approximation) tells us that the temperature from a given line-of-sight I in fact reflects the 
density, gravitational potential and velocities at the position xi = {rjo — ijr)1 and time rjji, as well as any time- varying 
gravitational potentials along that light-cone [27_ : 

e(x = 0,77o;f)= 0o + $ + ?-VH (xi,r7fl)-f / d7^{<^' ^ ^'){xum) , (4) 



where Vb is the baryon velocity potential. This equation shows that the primary (and ultimate) sources of anisotropics 
are the inhomogeneites in the matter and metric fields of the Universe (^Oj ^, ^ and Vb.) 

Now let's relax the assumption of instant recombination, but still require that the photons did not scatter again 
after decoupling. This is the same as saying that the visibility function is not assumed to be proportional to a delta- 
function 5(77 — rjn) anymore, but is still a positive, normalized function, highly peaked at the time of recombination. 
In that case the photons will carry the average temperature of the location where they last scattered, along with the 
local gravitational potential $ and baryon velocity Vb, and this signal will be affected by the time- varying metric 
perturbations only after that last scattering. Considering that the probability that a photon will scatter between 
some time 77' and the some time rj is given by the visibility function g{ri',ri), but the probability that they will not 
scatter anymore after the time 7;' is given by the total visibility ^{rj' ^rf), the line-of-sight solution becomes: 

e(i)(0,7;;f)=y" dij' [g[i^\T^)[eo + <^ + i-VVb]{xi,r,')+^{i^',ij){<^' + -^'){xuij')] , (5) 

where now xi = Arjl, with Ar/ = (77 — r/'). The superscript 1 is used to denote that this contribution is linear with in 
the total visibility 7 (as well as the visibility function, g = dj/drj). 

Equation ([5| tells us that, in a first approximation, to obtain the temperature anisotropics one should simply 
average the sources over the PLC xi, with weights given either by the visibility function (for the SW and Doppler 
terms) or by the total visibility (for the ISW term.) This approximation would correspond to the "one scattering" 
diagram of Fig. 1. 

The next level of complexity leads to polarization. Let's assume that the approximation above is still true for the 
temperature, but that photons can scatter a second time after decoupling. Then, the incident radiation at the location 
of that scattering will in general be anisotropic, simply because of the inhomogeneities in the Universe at the time of 
decoupling. If that incident radiation has a quadrupole, then Thompson scattering will excite the linear polarization 
degrees of freedom of the Stokes parameters Q and U. Clearly, then, polarization is of at least second order in the 
visibility, since it enters once when the photons first decouple from the matter, and then a second time when the 
anisotropic ensemble of photons scatter off free electrons, generating the polarization. 

The assumptions above are too simplistic, of course: as we go back in time the number of scatterings per Hubble 
time rise steeply, which means that the problem that must be solved is one of sucessive scatterings of a polarized, 
inhomogeneous and anisotropic temperature distribution which is, moreover, coupled to baryons and dark matter. 



The result of taking into account the anisotropy and polarization of the incident radiation in Thompson scattering 
leads to corrections to Eq. (l5| and to the generation of polarization [HI UHl HH HH HSJ |3H [33 • The polarization at 
any given point in space and time is best given in terms of the (dimensionless) spin +2 eigenstate combination: 



Q + iU 
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dv'g{ii',ii)\e2{x',7j')-V6p2ix',v') 



(6) 



where 62 and p2 are the quadrupole of the temperature and of the polarization, which we will define in more detail 
below. To these equations we should naturally add the perturbed Einstein equations, as well as the continuity and 
Euler equations for baryons, dark matter and neutrinos - see, e.g., [Sl ITTlfTB] . 

It should be evident from the symmetries of the problem that it is natural to break this system of equations into 
spherical coordinates with respect to the lines-of-sight I around an observer at the origin. For temperature, which is 
a scalar under rotations, the spherical harmonic decomposition reads: 



8(f , 77; i) = ^ Ohnix, r])Yhn{i) 



(7) 



Polarization, on the other hand, is such that the Stokes parameters Q and U change sign if we perform a rotation of 
TT over the line-of-sight, which means that they are components of a spin-2 field. In fact, the complex combination in 
Eq. (IgI was chosen such that it is a spin +2 eigenstate. Hence, polarization in this form can be expanded in terms of 
the spin +2 eigenfunctions, or spin +2 spherical harmonics |22lE5j : 



P{x,r];i) = ^ Pemix,ri)2Yirn{i) 

l>2,m 



(8) 



In Fourier space the dependence on I can be easily isolated, since e**^'^ = e^'^'^e*^'''^'. We then employ Rayleigh's 

expansion, e*'''^ = 4'K^^^i^j(^{kx)Yl^{k)Yirn{x), and the line-of-sight integrals determining anisotropics can be 
written in the form 1231 ESl : 



Qtm{x,rj) = 47ri' 



Ptm{x,-n) = 47ri^ 



(27r)3/2 

d^k 

(27r)3/2 



Y;,„(fc)0,(/c,77) 



'=-->;,„(fc)p,(fc,77), 



(9) 
(10) 



where the temperature and polarization momenta, 6^ and pi, are derived from the the geodesic equation for photons 
in the presence of Thompson scattering: 



ei{k,vi) = 9','\k,r]) + ^J\r,'gir]',r,)[92ik,r,')-V6p2ik,r,') 
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eo{k,fj') + <i>{k,7j') + Vb{k,v')f- 

dri' g{ri',ri) 02ik,ri') - V(ip2{k,r]') 



a(fcA77)2 
+ 7(ry', /;)($' + *')(&, Tj') \ j,{kAr,) , (12) 



[k^rif 



(13) 



For the derivation of the polarization term for the temperature quadrupole, including the radial function, see also 

^- . . 

Therefore, in this form the structure of the interactions is more clear than in the hierarchy of Boltzmann equations, 

but at the price of stating the problem in terms of integro-differential equations. This complexity is just apparent, 

though, since all higher-order source terms have angular momenta f < 2, and if we solve for the low ones, all the 

higher multipoles can be computed with the help of the integrals above |55] . 



III. CMB IN POSITION SPACE AND CAUSALITY 



The integral equations ( |11[ )-( |13| should be solved, in principle, for all momenta k so that the temperature and 
polarization anisotropics in Eqs. (|9I)-(10) can be computed. This includes the modes with k « Hq, which correspond 



to arbitrarily large wavelengths and can contribute to the zero mode of the fluctuations. Now, does this mean that 



perturbations outside the horizon can contribute anything to the CMB that we observer? In fact they don't: we will 
show that the temperature and polarization of the CMB which are observed at any spacetinie point (af, 77) only include 
information from inside the PLC of that point. This statement is true at each order in the visibility 7, and for each 
spherical mode (£, m) . 

The first step to recover the causal structure which underlies the temperature and polarization is to go from Fourier 
space back to position space. The most direct way to go back to position space without relinquishing the spherical 
harmonic decomposition is to use the fact that the Fourier- and position-space harmonics are simply related by a 
Hankel transform: 

/(2?) = E/iA^(^)^iA^(*) ' /W = E/iAf(fc)i"LM(fc), (14) 

LM LM 

fLM{x) = ^li'' J dkk^jLikx)fLMik) , fLM{k)^\jh~i)'' j dx x'' J L{kx) f lm{x) . (15) 

If /(if) is a real function, then the harmonic coefficients obey the conjugation relations f^„^{r) = (— l)™/f,-m(?') hi 
position space and fi^{r) = {—lY^™fi-m{r) in Fourier space. 

The relations above between fLAiix) and fLAiik) are quite remarkable: they tell us that the spherical harmonic 
phases do not mix at all. This is a consequence [39 of the fact that angular momentum is the same operator in 
position and in Fourier space, h = ix x d^ = ik x dk ■ It is worth noting that apparently this technique were first 
used in Cosmology in connection with redshift space distortions - see, e.g., PSI - HT] . In connection with the CMB, 
the spherical decomposition has been used in simulations [ 29| [5U] , and as a tool to study polarization from clusters of 
galaxies by [38] P]. 

A. 7^ term: Q^^' in position space 

Now we can easily substitute the sources in terms of this spherical harmonic decomposition into 6*^ . In the following 



subsections we show how this prescription can be extended to the remaining terms in the expressions (111 and (13). 
Let's then express the monopole ^O) newtonian potential <& and baryon velocity potential VJ, in terms of spherical 
harmonics Yl]\j(x), and use them to compute the temperature anisotropics at our location (assuming that we occupy 
the origin of the spherical coordinate system, at a? = 0), to first order in the visibility. Substituting the spherical 
harmonic decomposition in x into Eq. ( 12 1 for 9'^^\ inserting that expression in Eq. njj and integrating over d^k 



(which makes L = i and M = m), leads to: 



©12(0.^) = -/ dkk^ dr,' dxx^Si^{x,r,,7j')jiikAr,)jiikx), (16) 
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^0 Jo 



where we have collected the sources in the term: 

d 



Sijn{x,r],f]') = g{r]',T]) 



'0An{x, V') + ^emix, -q') - VbJmix, T]' ] , 



-j{v',f^){<^' + ^')e^{x,7i). (17) 



Notice that we have defined the sources in an unusual way, including a factor of the visibility function for the Sachs- 
Wolfe and Doppler terms, and a factor of the total visibility for the integrated Sachs- Wolfe term. The matter and 
metric fields are actually just functions of {x,ri'), but we include the extra dependence on 77 that comes from the 
visibility into the definition of the source term in order to simplify the notations. 



Now the integral over k in Eq. (16) can be computed, and in fact that happens to be exactly the orthogonality 



condition for spherical Bessel functions: 



f 

Jo 



vr 



dk k^ ji(kx) je{kx') = -x-^6{x-x') . (18) 
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This implies that Eq. (16) can be simplified to 

rn 



©£2(0,^) = / dTj' I dxS,^{x,r^,r]')S{x-Af]) (19) 



which is just the harmonic decomposition of Eq. ([5|. Notice that in Eq. (|5| x = /, so the gradient in the Doppler 
term can be written as a time derivative, which after integration by parts with the derivative of the delta-function 
becomes the derivative of the baryon velocity in the expression above. This warm-up exercise is useful to check that 
all sources which contribute to temperature anisotropics at this level come from the light-cone (its surface, in this 
case), which here appears explicitly as 5{x — Ai]). 

Another remarkable fact, which already shows up in this lowest-order approximation but which is true to all orders, 
is that the phases {i, m) of the CMB observables are the same as the phases of the sources. The only conditions for 
this to hold are, first, that we keep to linear perturbation theory, and second, that the optical depth is a function of 
time only. We will see next that this holds true to higher orders in the total visibility 7. 

B. 7^ terms in position space 



The expressions (11 l-(13l are integral equations for the temperature and polarization momenta 6e and pi. We can 
iterate these equations and organize the series into powers of the total visibility, similarly to what is done for the 
Dyson series of Quantum Electrodynamics |53] - except that the fields in the integral equations for the CMB are 
coupled to a set of ordinary differential equations (the Einstein, continuity and Euler equations for metric and matter 
perturbations.) 

In the previous subsection we computed the first term of this series, which is of order 7^ (since the source term is 
itself linear in 7). Consider now the next terms of this series, which are of order 7^. 

For polarization, we have that, to order 7^, the only term which contributes is the temperature quadrupole to 



order 7^. By substituting Eq. (12) into Eq. (13 1 and expressing the sources in terms of their spherical harmonic 



decompositions, like was done in the previous subsection, we obtain after some algebra: 






Pi'JiO^v) = -^J^^^jy^^'div',!)/^' d^" J^ dxSU^,v',ri") (20) 



2 • /, ^ JiJkAr]) 
(fcA77)2 



X / dkikxfMkx)^-l^^j2ikAri'), 



where Arj' = rj' — rj" . A simplified version of Eq. (20) was first written in this form by [38', for the case of the 
polarization from a galaxy cluster - where the visibility function after decoupling can be thought of being proportional 
to a Dirac (5-function, gc{rj,x) — fic5{x — hArj), where fic is a cluster's optical depth and n is the line-of-sight to the 
cluster. 

Physically, this contribution to polarization corresponds to the photons that decoupled at some radius x and time 
77", are scattered at time 77', and then end up as a polarized beam at time 77. The geometry is shown in Fig. 2. We will 



show now that the integral over k at the end of Eq. (20) has exactly that meaning: it vanishes unless the distances 
X, Arj and Arj' form a triangle, and such a triangle does not exist unless the sources are inside the PLC of the last 
scattering point, and the point of last scattering lies inside the PLC of the observation point. We will give a general 
expression for integrals such as this in Appendix El but here is the result [28| : 

„2 /-oo 

W!{n,r2;r3) = ^ / dk u(kri) ],{kr2) J2{kn) (21) 



where: 



- — \Pe (cosQ;i2)sin^ai2 , 
4 r2rt 



cosai2^ ^ ''' ' (22) 

2rir2 

is the cosine of the angle between the sides ri and r2 in the triangle of sides ri , r2 and r^ . The window function Wf is 
zero if that triangle does not exist, which means that it is nonzero only if the following set of conditions are satisfied: 

ri <r2+r3 , r2 < rs + ri , r^ < ri + r2 , (23) 

or, equivalently, \ri — r2\ < r^ < ri + r2 - see the left panel of Fig. p] This window function was first computed in 
Ref. [38]. We have plotted some cuts of those window functions in the Appendix, Fig ??. 



Hence the momentum integral in Eq. (25) is W^{x, Ary, Arj'), which means that it is a spacetime window function 



that vanishes unless the inequalities above are satisfied. What this result implies to our Eq. (20) is that a source at 
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Figure 2: Left panel: if the triangle of sides ri, r^ and rs exists, then the window function Wf is non-zero. Middle and right 
panels: spacetime diagrams for the scatterings, where time runs up, and space (radial coordinates) runs horizontally from the 
center to the sides. The diagrams show sources which decoupled at r{\ then scatter off free electrons at time 77', and end up 
contributing to the polarization P\J^ at time (r\). The minimal and majcimal values of x for which the sources can contribute 
to polarization are given by Xmin = | A77 — A?;'] and Xmax = A77 + Ar;' = 77 — 77". 

radius x and time r\' can only contribute to PfJ^ (77) if | A?7 — A77'| < a; < A77 + A77' = rj — rj" . As the right panel of Fig. 
2 shows, this spacetime window function limits the contributions of the sources to the PLCs of the last scatterings, 

and the scatterings themselves to the PLC of the observation point. Obviously, this means that the sources that 

(2) 
contribute to P^^ must also be inside the PLC of the observation point - which in this case is 77 — 77". 

Let us summarize this result for the lowest-order contribution to polarization in position space: 



p^Mv) = 




JO 



dr]" 



n-v 



dxSim{x,v',ri")W!ix,ATj;ATj') 



(24) 



Notice that, among other restrictions, the window function above cuts off the spatial integral at radius x at rj — rj" , 
which is where the PLC of the observation point lies at time 77". 

In ^38] Eq. ( |24[ ) was used as the starting point to show how to invert the polarization data from galaxy clusters in 
order to reconstruct the three-dimensional map of the sources at the time of decoupling. That constitutes the solution 
(to lowest order in visibility) to a conjecture by Kamionkowsky and Loeb about how to get around cosmic variance 
using cluster polarization [43J. 

Therefore, just as happened in the previous subsection with &£„[, the integral over k in effect guarantees that the 
sources will only be taken into account if the scattering processes happen on the light cones. As before, the integral 
over X in Eq. (20 1 is cut-off at Xmax = V ~ v" ^ eliminating sources which lie outside the PLC of the observation point 
(0,77). But now there is another feature: since Xmi7i = IA77 — A77'|, sources which were too close to the observation 
point at time 77" also cannot contribute to the CMB at time 77 if last scattering happened at time 77' - see the middle 
and right panels of Fig. [2J This additional constraint on the volume of the PLC which is integrated applies for times 
of last scattering (77') which are both close and far from the observation time 77. 

Before we turn to the order 7"^ terms, we write down the order 7^ contribution to the temperature anisotropics, 
which comes from inserting 62 into Eq. ( |ll| . After a calculation very similar to the one done above for polarization, 
we obtain that: 



0f2(o.^) 
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(25) 
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The integral over k on the last line of the previous equation can be recast in terms of W^ if we use the recursion 
relations for the derivatives of spherical Bessel functions: 



.Te_ 



Z'-— [z=^^jf(z)] =±J£zpl(z) 



dz 



(26) 






Figure 3: Diagrammatic representation of the series in the visibihty. In all diagrams, the filled dot corresponds to the point 
and time of observation of the CMB temperature and polarization. The 7^ term (the two lines beginning and ending at the 
same points) corresponds to a window function 5(ri — r2). The j'^ a nd 7 '^ terms (respectively the triangle and 4-side polygon) 
correspond to the window functions W'i of Eq. ( 21 \ and Wl of Eq. ( 28 1. For all diagrams the window functions vanish unless 



the sides are such that the polygon can be closed. The main physical implication is that sources (which by default are located at 
n with respect to the observation point) outside the PLCs of the last scatterings do not contribute to the CMB. In particular, 
sources outside the PLC of the point of observation are thrown out of the integration of the physical observables. See the 
Appendix for a full discussion of these geometrical properties. 



The niomcntum integral then becomes 



{xAriY 



dx dArj 
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which shows that, just as happened for polarization, the contribution of order 7^ to the temperature is also modulated 
by the same spacetime window function. 

C. 7'^ and higher-order terms in position space 

The iteration process is trivial, but the higher-order terms become lengthy. We present the result for the simplest 
7'^ term that contributes to polarization: 



-P/m(0'^) = -^^\\to — ?aT / drj' g{ri',ri) / dr]" g{r]" ,r)') / drj'" / dx Si„,{x,t]" ,r]"') 
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,0 (fcA77)2 {kAj^'Y J2ikAr,") , 

where A77" — ry" — ry'". We will show in Appendix ^ that the integral over k can be resolved into: 



(27) 
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2 J r^ 4 r2r-^ 



(28) 



where now r is defined as the common side of two triangles, of sides (?'i,r2,r) and (r3,r'4,r), such that r^ = rf+r2 — 
2rir2 cos ai2 — r^+ r1 — 2r^r4^ cos a34 - see Fig. p] It is trivial to show that all the remaining terms of order 7"^ that 
contribute to the CMB temperature and polarization can be written in terms of the window function Wf by using 
the recursion relations of spherical Bessel functions. 

The window function vanishes unless all four sides can form a (flat) polygon, so inequalities similar to those found 
for wf apply here as well: 



'rs<fi + 'ri+ r2 



r2<r'i + ri + ri , 



(29) 
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Figure 4: Structure of the successive PLCs for two (left diagrams) and three (right diagrams) scatterings. The values of x for 
which the window functions Wi{x,'ri;ri') and Wf{x^r\\r\ ,vl'^ are non-zero are indicated as the thick lines at the bases of the 
cones. Here the radius x should be measured from the origin of the spatial coordinates, which in these spacetime diagrams lie 
along the dashed lines. Notice that the region in the vicinity of the PLC of the observation point (the outermost PLC) always 
contributes to the integral, the origin is always excluded, and intermediate regions may or may not contribute to the CMB 
observables. 

These conditions imply once again that the sources can only contribute to the observable at time r\ if each successive 
PLC lies inside the previous PLC, all the way from the observation point back to the sources. The constraints on 
the positions of the sources are now more complicated than the case of two scatterings, but the diagrams in Fig. |4] 
show which values of x (measured from the origin, at the middle of the bases of the cones) are allowed by the window 
functions. In three spatial dimensions, these ranges of radii correspond to concentric spherical shells - in the case 
of two scatterings, there is one spherical shell; in the case of three scatterings, two spherical shells; and so on. The 
outermost spherical shell always includes the edge of the PLC of the observation point. 

It is easy to recover Wf from Wf by taking ra -^ 0. In fact, we can also recover the orthogonality condition for 
the spherical Bessel functions, Eq. (18 1, from Wf by taking r-j, — )■ 0. This can be made by noticing that, in Wf the 
limit r3 -> leads to r-i = ri, and in W/ the limit r3 — >■ leads to r4 = r. Then, using the expansion of the Bessel 
functions for small arguments, J2{z) ~ z^/15 + ©(z^), we obtain that: 



lim VKf(ri,r2;r3) = 

ra-fO 



15 
1 



r^ ^(5(ri - r2) 



lim Wi{ri,r2;r3,ri) = — x W^{ri,r2;r3) . 

ra—yO io 



(30) 
(31) 



The first identity just shows that the PLC delta-function is the spacetime window function for two scatterings (order 
7^.) The second identity can also be verified by noticing that in the ra — > limit, cosa34 is unscontrained, so we can 
integrate out the dependence on that angle, which gives a factor of 2/15. 

The procedures outlined above can be extended once again to the next order in the visibility, 7**. The geometrical 
interpretation is given by the rightmost diagram in Fig. [3] Again, we see the role played by the spacetime window 
functions, of regulating the volume of the PLC in each scattering so that the information from the sources is propagated 
causally all the way to the observer. Notice also that by taking ra — >■ we recover W* - in fact, all the spacetime 
window functions obey the relation: 



lim W"{ri,r2;r3, 

1-3 ->-0 



,rn) 
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D. Interpretation of the series in 7 

In this Section we have shown, first, that the contribution of order 7^ from sources at positions x and times of last 
scattering 77' to the temperature which is observed at time t] is modulated by a delta-function, 5{x — A77) - and this is 
nothing but the surface of the PLC of the observation point, x = Ary = rj — rj'. Then we showed that the contributions 
of order 7^ to temperature and polarization are modulated by the spacetime window function W^{x, Arj; Ar/), where 
now At] denotes the interval between observation and the last scattering, and Arj' denotes the interval between that 
last scattering and the second-to-last scattering. We presented the simplest contribution of order 7^, which sums 
over sources modulated by the window function W^, and indicated how the order 7'* term is also modulated by a 
window function W^. Notice that we have been labeling each term of the series in terms of 7, which in this notation 
stands both for the total visibility 7 = e~^ and for the visibility function g — d'y/drj. However, even though the total 
visibility is always smaller than unity, the visibility function is highly peaked at the time of recombination and during 
reionization, so it would be incorrect to characterize our series over visibility as a perturbative series - it is rather 
more like an asymptotic series. 

The series over visibility has the diagrammatic representation shown in Fig. |3] Perhaps it is not so surprising that 
the spherical Bessel functions play a key role in regulating the volume inside the PLC, since they are associated with 
the matrix elements of E3 ~ the Euclidean group in three dimensions [H]. E3 is a non-compact group consisting of 
the set of transformations under which distances are invariant, which may explain why the spherical Bessel functions 
modulate the (invariant) distances on each equal-time hypersurface. 

Another check on our results is the fact that all the window functions are real. Since the sources S{x, rj, rj') are 
themselves real, this means in particular that P^^ = (— l)™P£__m - i-e., the polarization that is generated from 
scalar perturbations is made up purely of _E-modes. If we had included gravity waves or lensing, these window 
functions would have acquired an imaginary piece as well, which would have "magnetic" (instead of "electric") parity 
and therefore generate S-modes. 

The spacetime window functions to any order in 7 vanish unless each scattering lies in the PLC of the following scat- 
tering, all the way from the sources to the observation point. The order 7^ contributions to CMB observables hold the 
information from the iV-th last scatterings, and are modulated by the window function W^'^^(x, Arj] Arji, . . . , At^tv-i)- 
In terms of the diagrams of Fig. [3] there are always two spherical Bessel functions of order £, corresponding to the 
source (x) and the last scattering before observation (Ar/), and A^ — 1 spherical Bessel functions of order 2, corre- 
sponding to the iV — 1 intermediate scatterings for which the quadrupoles of the temperature and polarization served 
as sources in the iterations of the line-of-sight integrals. The computation of these spacetime window functions, as 
well as the restrictions they impose on the sources that can contribute to the CMB (in particular the fact that they 
all vanish outside the PLC of the observation point) are shown in Appendix [A| 

Therefore, in position space the line-of-sight integral equations have a simple interpretation in terms of successive 
Thompson scatterings happening along the successive PLCs. These expressions should still be coupled to the Einstein, 
continuity and Euler equations through the temperature dipole and quadrupole, but those are local equations so their 
causal structure is trivial (nevertheless, the Green's function in position space for the cosmological matter and metric 
perturbations can also reveal very interesting features |34[[35|.) 

Going back to the hierarchy of Boltzmann equations, which are obtained directly from the line-of sight integrals 
in Fourier space, there is one full hierarchy which holds separately for each mode {k,i,m}. Hence, what we have 
shown is an explicit demonstration that the line-of-sight formalism can also be seen as the solution to an initial values 
problem, where the initial conditions only need to be specified inside the PLC of the observer, at some initial time 
when the visibility was small enough that the series converges quickly. In other words: we have explicitly shown that 
the line-of-sight integrals in Fourier space are essentially the Fourier transform of a retarded Green's function for the 
CMB observables. These retarded Green function for anisotropics in position space are expressed at each order in the 
series over visibility through spacetime window functions, which are ultimately the objects responsible for enforcing 
causality and regulating which inhomogeneities are able to affect anisotropics, and how. 

The results of this Section have an interesting connection with schemes to simulate constrained maps of CMB 
temperature and polarization J29H33| . These simulations are extremely important to predict the types and levels of 
non-gaussianity in the CMB which are generated at the very early Universe, since distinct inflationary models can 
be differentiated on that basis |45| . In those simulations, the temperature and polarization transfer functions, which 
transform the inhomogeneities (the sources) into anisotropics, are computed numerically assuming some visibility 
function. Moreover, some approximations are typically made, such as considering only the curvature perturbation in 
the source term. Our expressions, on the other hand, are exact and analytical, but it is not immediately clear how (or 
if) they could be used to facilitate a simulation. However, our results show that the transfer functions of 1291150] are 
made up of invariant (or geometrical) pieces which can be factored from the purely time-dependent visibility function. 
These invariant parts are given by our spacetime window functions, which do not depend on the cosmological scenario 
or on the history of recombination. The transfer functions relevant for the simulations can, therefore, be obtained by 
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integrating our spacetime window functions over time with some visibility function. 

IV. CMB WITH THE FOURIER-BESSEL EXPANSION 

We have seen that the Hne-of-sight integrals for the CMB temperature and polarization, when framed in position 
space, lead to spacetime window functions that constrain the positions of the sources S{x) that are eventually inte- 
grated over. These constraints, valid at each order in the visibility, imply in particular that the sources must all lie 
inside the PLC of the observer that measures the CMB. But this is not the only constraint that is relevant for this 
problem. 

The physics of recombination is such that the visibility is exponentially small for z ^ lO'^. Since all the terms of 
the series over visibility are linear on the sources, the fact that they are always multiplied by powers of the visibility 
implies that the sources which are relevant in the line-of-sight integrals are those that lie inside the PLC of the observer 
at some time late enough that the visibility is non-vanishing. So, the sources which are actually summed into the 
CMB anisotropics are not simply the sources in the PLC of the observer, but the sources in the PLC of the observer at 
times such that the visibility is not totally negligible. Since for a typical flat ACDM cosmology the lookback distance 
from today to those early times saturates near R ~ 5Hq^ for z ^ 10^, it makes little difference whether we choose 
that initial instant (when 7 — > and g — t' 0) to be z = 10^, 10^ or 10^". 

This means that the integration of the sources is not simply limited to the PLC of the observer, as implied by the 
window functions Wf^ , but that the relevant spacetime volume is that of the PLC, cut-off at some initial time rji 
such that 7(77, rji) and g{ri, rji) are sufficiently small. Since the unperturbed spacetime is symmetric around the source, 
there is in effect a maximal radius away from the observer, R = rj — rji, such that outside that radius, the sources are 
effectively zero by virtue of the powers of 7 that multiply them at each order in the series. Notice that we could even 
choose R to be bigger than this lookback time, but that would unnecessarily include sources which are eventually 
discarded in the integration of the physical observables. 

Hence, we can fix some boundary R and set all fields to zero at that boundary and beyond: the result of solving 
the initial values problem through the line-of-sight integrals will be exactly the same, order by order in the series over 
the visibility. 

If that is the case, then we should ask what would be the best way to represent the sources, considering that they are 
zero at and above some radius R from the origin. The appropriate expansion in that case is clearly the Fourier-Bessel 
series [25], for which the fields are expanded in spherical harmonics and spherical Bessel functions, but instead of the 



continuum of momenta that appears in Fourier space, Eq. ( 15 1, the modes are discretized: they are given by the roots 



of the spherical Bessel functions. A function that obeys Dirichlet boundary conditions at r = i? is expanded as: 

CO 

J{x) = ^ ftm{x) Yi>r,i{x) = X! X! ■^*^" kihfX) Yim{x) , (32) 

im £m i—1 

where the last sum is over the i-th root of jg, so that for all i's the Bessel functions vanish at the boundary, jeikuR) — 0. 
The coefficients fum can be obtained by using the orthogonality relation of the Fourier-Bessel basis: 

f^ 1 

/ dzz'^ji{quz)ji{qjiz)^-[je+i{qu)f6ij, (33) 

where qu is the i-th root of ji{z). This expression leads inmiediately to: 



2R-3 j-R 

Ji+iiheR) Jo 

where k^ — qu/R. 

The most important feature of the Fourier-Bessel series is that its basis functions, like the plane waves of the Fourier 
expansion, are eigenvectors of the Laplacian operator in flat space, 'S/'^ ji{kiix)Yim{x) = — kfgji{kux)Yimix). This 
means that the Einstein, continuity and Euler equations for matter are exactly the same as in the usual Fourier 
expansion - except that the momenta ku are now discretized. Moereover, the lowest eigenmode of the Fourier-Bessel 
series that can contribute for a given multipole £ is ku ^ £/R. 

A. Fourier-Bessel modes of CMB observables 

Now we can go back to the CMB and translate the equations and methods presented in the previous Sections to the 
Fourier-Bessel basis. The most important difference between the two expansions can be grasped by comparing Eqs. 
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(15 1; (32 1 and (34 1. Basically, when going from the Fourier basis to the Fourier-Bessel basis, the angular dependence 



is still expressed in terms of spherical harmonics, but the radial coordinate is expressed by a sum, not an integral 



dkk' 



R-'H^ 



(35) 



which just tells us how to go from the continuum of momenta appropriate for fields in M"^ to the discrete tower of 
momenta kn that encapsulates all the information for fields limited to the finite volume inside a sphere of radius R. 
In particular, this means that now the CMB temperature and polarization are given in terms of the momenta not 
by Eqs. ([9]) - ([lO|, but by: 



■i 
Plmiv) = ^Piemiv) ■ 



(36) 

(37) 



In these sums Oum = (^i,Um ^^nd Pum = Pi.Um £^re given by the solutions of integral equations analogous to Eqs. 
(11 1-(13), which for each multipole L and for each mode {i£m} read: 
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(fc^^A7,)2 



jdhiAvi) , (38) 
(39) 

(40) 



and where the sources in Eq. (39), which were defined in Eq. (17 1, have been expanded in the Fourier-Bessel basis 



as well. Notice that only the generalized modes {L;i,£,m} with L <2 really need to be calculated from the integral 
equations (or, equivalently, from the associated Boltzmann equations), since all the higher modes (L > 3) can be 
computed from the former, and only the pieces L — i actually get summed into the CMB observables, Eqs. (36l-(37l. 



This is just a restatement of the fact that the sources of anisotropics are the matter and metric inhomogeneities, 
plus the dipole and quadrupole of temperature and polarization - which, again, underpins the vast superiority of the 
line-of-sight formalism compared to the full hierarchy of Boltzmann equations. 

An important check of consistency is to reobtain the temperature and polarization anisotropics in position space 
that were derived in Section III. To first order in 7, the temperature anisotropics are given by: 



Qilliv) = E^- 



(1) 



(41) 



drj' I dx X 



2R 



-3 






where the source term Sim{x,r],r]') was defined in Eq. ( |17[ ). But now the infinite sum in Eq. (41 1 can be resolved 
through the use of the orthogonality of the Fourier-Bessel basis in target space (actually, this is a completeness relation 
- see [28j, Cap. XVIII), which is just the Fourier-Bessel counterpart of Eq. ( 18 1: 



E 



Mikun)ji{kar2) _ 1 3 2 ., ^ ^ 

j/+l(fe«rt) 2 



(42) 



This identity leads then automatically to Eq. (19 1, which shows that to order 7^ the Fourier and Fourier-Bessel 
descriptions are identical. 

To order 7^ it is less obvious that one obtains the same anisotropics as related to the position-space fields, but it 
is true nevertheless. To see that, take the simplest term - the contribution to polarization that comes from the order 
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7^ quadrupole: 

Pi^iv) = ^PiMrn (43) 




2V /"' /"' 



drj' g{ri' ,7]) / dr;" / dx 



;a: 



Notice that now, making ri = x, r2 = A?7 and r^ = Arj' we have on the right hand side a spacetinie window function 
given by: 

Although we haven't been able to prove mathematically that this Fourier-Bessel window function is indentical to the 



Fourier window function of Eq. (21 1, we have checked numerically that they are identical - including the factor of tt 



which relates the phase spaces of the two basis functions. Hence, the lowest-order contribution to CMB polarization 



that results from using the Fourier-Bessel representation is again given, precisely, by Eq. ( 24 ) . 

In fact, we can prove (see the Appendix) that the not only the two window functions W^ and W^ are equal, but that 
all the window functions derived in Section III are identical to the window functions that arise in the Fourier-Bessel 



expansion, if a certain generalization of the orthogonality relation, Eq. (42 1, is valid 



ji'{kiiri)ji,{kiir2) _ 1 



^ Ji+iiheR) 2 



= 7;R r^ Sin-r2). (45) 



i=l 



We have checked numerically that this relation seems to hold true for a range of i and £' , and for any arguments 
< 7-12 < 1, but as far as we know this has not been proven anywhere in the literature about Bessel functions - even 
though it clearly is a fundamental tool for relating quantities in the Fourier and in the Fourier-Bessel expansions. 



We have also checked that the integral equations (38l-(40) lead to the usual hierarchy of Boltzmann equations 



[161 1231 E5] . where now there is one full independent hierarchy for each mode {i£m}. It is curious that, while in the 
usual Fourier analysis what generates the hierarchy of Boltzmann equations are the recursion relations of the Legendre 
polynomials in the angular dependence k ■ I, in our case the generators of the hierarchy are the recursion relations of 
the radial modes - the spherical Bessel functions. Since the two special functions are intimately related by Rayleigh's 
expansion of the plane wave (which is ultimately what regulates the line-of-sight integrals) , it is indeed natural that 
both basis could be used to generate that hierarchy. 

The crucial difference between the Fourier-Bessel series and the usual Fourier analysis lies in the discrete momenta 
kii that can contribute to the observables Qim and Pi^ in the Fourier-Bessel expansion. Critically, in this discretized 
series the first mode to contribute at each multipole is ku ~ £/R. Another feature is that the total number of modes 
which one needs to compute to obtain anisotropics up to some Imax is Nmax ~ ^max/^ (fo'^ £max ^ 10). 

There is, however, an apparent drawback of the Fourier-Bessel basis: although it clearly is a superior method to 
solve our sort of initial values problem when compared to the Fourier basis, when the underlying spatial fields are 
Gaussian the coefficients of the Fourier-Bessel series do not obey simple statistics like those of the Fourier modes. We 
will now turn to these issues. 

B. Statistics and power spectra in the Fourier-Bessel basis 

Take a Gaussian field f{x) in M.^, which is expanded into spherical harmonics in position and in Fourier space as in 



Eqs. (14). If homogeneity and isotropy are unbroken, the two-point correlation functions in position and in Fourier 



space can be expressed in terms of the radial functions as: 

{fe^ir)fi'ra' ir')) = ^«' S„,,n' ^/(r, /) , (46) 

{fim{k)f;,^,{k')) = Sw S,n,n' fc"' Pf{k) 6{k - k' ) . (47) 

Here the advantage of Fourier space becomes evident: translational invariance implies that the covariance matrix of 
the Fourier modes is completely diagonal. The relationships between the position-space two-point correlation function 
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and the power spectrum are given by: 

ilir,r') = - / dkk^Mkr)jeikr')Pf{k) , (48) 



and, conversely, by: 



2Je(Mlcf/ 



Pf(k)^J^ d"^"- J^ft^'-y)- (49) 

From the identity above it is also evident that the two-point correlation function has a lot of redundant information, 
since many different traces of it can lead to the power spectrum. If the field / is Gaussian, these correlation functions 
are the only non-trivial statistical momenta of that field's distribution function. 

We want to obtain the corresponding relations for the modes of the Fourier-Bessel series. This can be easily achieved 
by taking fields f{x) in M"^ and passing them through a radial window function W{r), such that W{r > R) — - e.g., 
the tophat window function: 

W^''ir) = eiR-r) , W^"ik) = ^^nikR), (50) 

where 0{x) is the step (Heaviside) function. In this way, the filtered functions will obey the boundary conditions, 
fiR) = WiR)fiR)=0. 

In Fourier space, the effect of a window function is to couple the different modes: 

/(fc) - / J^ W{k q) m . (51) 

Hence, the Fourier transform of the filtered function acquires non-diagonal correlations /c-space: 

imr (k')) =I^Wik~q) W{k' - q) Pf{q) . (52) 



For a purely radial window function W{r), we obtain with the help of Eqs. (15 1 that the spherical harmonic 
components of the filtered function are: 



fimik)= I dq q' /,„ (q) Wiik,q) , (53) 

where We{k,q) is the (symmetric) mode-coupling kernel of the radial window function: 

2 f°° 
Wi{k,q)^- drr^W{r)jfXkr)Mqr). (54) 



Hence, it is clear that if W(r) -^ 1 then Wi{k,q) — > q^^5{k — q) and we recover the Fourier modes of the M.^ field. 
For a generic radial window function, however, there will be mixing of modes, and the covariance matrix will be 
non-diagonal. The filtered spectrum is then related to the physical power spectrum through: 

{fimik)fUk')} = r dqq^W,{k,q)W,{k',q)Pf{q). (55) 

Jo 

In particular, for the Fourier-Bessel modes, which are related to the (filtered) spherical modes in Fourier space by 

/i£„ = V^TT'i'^ Ji+iikiiR)femikii), this last identity implies that: 

{him f Jim) = -2 (I. m Al (I. m / dqq^We{hi,q)Wi{kji,q)Pf{q) 

J l+lVf^il -tl) Ji^l[K ji It) Jq 

^ 2 R-^ji(kur) r"^ ^^,^,2 2R-^u{kjir') ^ 2 

Je+i{kuR) Jo i}+i{kjiR) -^ Jo 

Ji+i{kiiR) Jo Ji+iikjiU) 



^^^ 72 ;rm / ^^'^ :2 ;; ^ X^ / dqq'Pfiq)Mqr)Mqr') 
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where we have used the tophat window function from the first to the second line. These expressions show how to 
compute the covariance of the Fourier-Bessel modes from either the power spectrum or from the two-point correlation 



function in position space. It is also useful to obtain the equivalent of Eq. (48 1 in the Fourier-Bessel representation. 
By the completeness relation, Eq. (|42l, it is easy to see that: 



i J 

where 6{r) is the step (Heaviside) function. 

The two-point correlation function (as opposed to the Fourier spectrum) is more directly related to the physical 
observables, since it remains invariant as long as we keep within the causally accessible region (i.e., r < R and r' < R.) 
However, Eq. ( [SG] ) also tells us that the two-point function of the Fourier-Bessel modes has non-diagonal terms (albeit 
only in fc-space.) This does not pose a problem, because the spectrum is not really an observable: it can only be 
estimated (with exactly the same tools and assumptions as usual) from observables such as the temperature and 
polarization maps, as well as their derived products such as the angular power spectra. And since we saw in the 
previous section that the observables retain exactly the same relations to the sources of anisotropics as they do in the 
usual Fourier expansion in K^ , we conclude that the Fourier-Bessel expansion fulfills all the requirements to faithfully 
express the physics of the CMB. 

C. Angular power spectra 

Most of the useful cosmological information that we get from the CMB comes from the angular power spectra, 
because of their simple relationshi p w ith the Fourier power spectrum. For a function /(a?), the angular spectrum at 
radius r can be defined from Eq. (46 1, as C}{r) — ^^{r,r). Below we show that in the Fourier-Bessel expansion the 



angular power spectra assume exactly the same values as they would if we did not assume that the space was limited 
to the sphere r < R. 

The argument is simplest for the temperature anisotropics to order 7^, and generalizes in a trivial manner to the 
higher-order terms. The angular power spectrum for temperature in the Fourier-Bessel case reads: 

(0l2©^'t) = (E ^'ill E ©i'^t') = Cr^'Hv) Su- 6^^, . (58) 



Using the orthogonality conditions of the Fourier-Bessel basis in target space, Eq. (42 1, we obtain: 

Cr^'Hv) = f dv[ f dv'2 {SU^ = A7^,,rj[)SUx = A^, V2)) (59) 

= {SemSern)PLC ■ 

But this is exactly the usual result: to lowest order, the temperature angular power spectrum is given by the average 
over the PLC of the (angular) two-point angular correlation function of the Sachs- Wolfe, Doppler and integrated 
Sachs- Wolfe source terms, properly weighted by the visibility. For the j'^ and higher-order terms, the procedure is 
precisely the same and leads back to the same relation between the angular power spectra and the sources as happened 
in position space, where the spacetime window functions regulate which sources contribute to the anisotropics in the 
space and time integrals. Hence, we have shown that not only the observables (the temperature anisotropics), but 
also that the statistics of the angular power spectra in the Fourier-Bessel expansion are identical to the usual case of 
the Fourier expansion. 

D. Fourier v. Fourier-Bessel 

The identity between the spacetime window functions at all orders in 7 implies that the source and the observables 
(the temperature and polarization maps, or equivalently their spherical harmonic components Qem ^^nd Pem) £^re 
related in exactly the same way in the Fourier-Bessel basis and in the Fourier basis. The statistics of the angular 
power spectra, therefore, are also related to the statistics of the underlying matter and metric fields in precisely the 
same way in the two representations. 

Hence, the Fourier-Bessel basis is, in some respects, completely equivalent to the Fourier basis: it represents the 
same physics and it expresses the same observables as its Fourier counterparts - it even has the same statistics. 
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However, in at least one respect the Fourier-Bessel basis is superior to the Fourier basis: it has a precise prescription 
for the discretized tower of modes that contribute for the observables at each niultipole. These modes take into account 
exactly the relevant pieces of information from the sources, the ones that propagate from the initial value surface to 
the physical observables - no more, no less. And the statistics of the power spectra, as we have demonstrated above, 
is related in precisely the same way to the statistics of the (presumably Gaussian) matter fields, just as happens in 
the usual analysis in Fourier space. 

Finally, as an initial-value formulation the Fourier-Bessel expansion is vastly superior to the Fourier representation 
because it does not waste any resources keeping track of irrelevant variables such as super-Hubble modes or modes 
which trace out of the observable. All the information is encoded in a discrete series of momenta, and we do not have 
to guess how to subdivide the Fourier space in sufficiently small pieces in order to sample the observables we want to 
compute - the Fourier-Bessel modes already provide the unique, optimal choice. 

V. CONCLUSIONS 

In this paper we have shown how causality constraints in position space regulate which sources of anisotropics (the 
matter and metric perturbations) can contribute to the CMB. This causal structure is manifested order by order in a 
series of terms corresponding to the number of interactions that photons experienced over the past light-cone of the 
observer - or, equivalently, a power series on the visibility 7 = e~''. 

When expressed in position space, the line-of-sight integrals acquire an intuitive interpretation in terms of scatterings 
over the light-cones of the successive scatterings, all the way from the sources to the point and time of observation. 
In particular, we find that, in position space, only the sources of anisotropy that are inside our past light-cone are 
taken into account. This statement is exact to all orders - as it should be, since the causal nature of the propagation 
of photons is the key ingredient in the line-of-sight integrals from which we started. 

At each order in the power series on the visibility, the sources are weighted by spacetime window functions. These 
window functions can be complicated for a high number of scatterings, but they all obey a very simple rule: they 
vanish indentically unless some extremely simple sets of inequalities are satisfied. These inequalities have a simple 
geometrical interpretation: if the position of the source and the radii of the light-cones of the interactions cannot form 
a flat polygon, the spacetime window functions vanish. An interesting question which we did not have time to address 
is at what number of scatterings prior to free streaming we can cut off this series so that the error in the temperature 
distribution is, say, of order 1%. 

One of the implications of these causality constraints is that, whatever the properties of the Universe outside a 
limiting radius R, the source fields do not propagate to the CMB observables - and the line-of-sight integrals both in 
Fourier space and in the Fourier-Bessel expansion retain this property. In practice, this means that we can use the 
Fourier-Bessel framework to compute the CMB - even though it puts the Universe in a "spherical box", and discards 
all the information outside of that box. 

In the Fourier-Bessel basis, the fields are decomposed in spherical harmonics and a series of discrete eingenmodes 
ki£ (as opposed to the continous modes of Fourier space.) The first eigenmode for each multipole £ is ku ^ i/R. CMB 
observables are exactly the same as in the Fourier basis - but the Fourier-Bessel basis is optimal, in the sense that it 
does not keep track of irrelevant modes, only the ones that contribute constructively to the physical observables. 

The previous discussion implies that our results and methods are suitable for analytical and numerical studies of 
CMB temperature and polarization maps in models with large-scale inhomogeneities, statistical anisotropy or non- 
gaussianities of any kind kind |2S]-(32]. Our results may be useful also in simulations of the CMB in the presence 
of topological defects. It is not clear whether the methods described here can be employed to study models with 
non-trivial topology [33], since in those cases the Fourier or Fourier-Bessel basis functions may not be eigenvectors of 
the Laplace-Beltrami operator. 

Acknowledgements 

The authors would like to thank J. C. A. Barata for many conversations on the properties of Bessel functions, and 
to Mathias Zaldarriaga for useful comments. This work was supported by FAPESP and CNPq. 



[1] G. F. Smoot ei ai, Astrophys. J. 396, LI (1992). 

[2] WMAP, C. L. Bennett et ai, Astrophys. J. Suppl. 148, 1 (2003), astro-ph/0302207. 



18 

[3] WMAP, D. N. Spergel et al, Astrophys. J. Suppl. 148, 175 (2003), astro-ph/0302209. 

[4] E. Komatsu et ai, (2010), 1001.4538. 

[5] WMAP, L. Page et al, Astrophys. J. Suppl. 170, 335 (2007), astro-ph/0603450. 

[6] QUaD, . M. L. Brown et al, Astrophys. J. 705, 978 (2009), 0906.1003. 

[7] H. C. Chiang et al, Astrophys. J. 711, 1123 (2010), 0906.1181. 

[8] S. Dodelson, Modem Cosmology (Academic Press, 2003). 

[9] V. Mukhanov, Physical Foundations of Cosmology (Cambridge University Press, 2005). 
[10] R. Durrer, The Cosmic Microwave Background (Cambridge University Press, 2008). 
[11] P. Peter and J. -P. Uzan, Primordial Cosmology (Oxford Univ. Press, 2009). 
]12] S. Chandrasekhar, Radiative Transfer (Dover, 1960). 
113] P. J. E. Peebles and J. T. Yu, Astrophys. J. 162, 815 (1970). 
114] J. R. Bond and G. Efstathiou, Mon. Not. Roy. Astron. Soc. 226, 655 (1987). 
115] J. R. Bond and G. Efstathiou, Astrophys. J. 285, L45 (1984). 
116] C.-P. Ma and E. Bertschinger, Astrophys. J. 455, 7 (1995), astro-ph/9506072. 
117] W. Hu and N. Sugiyama, Astrophys. J. 471, 542 (1996), astro-ph/9510117. 
118] W. Hu, N. Sugiyama, and J. Silk, Nature 386, 37 (1997), astro-ph/9604166. 
119] U. Seljak and M. Zaldarriaga, Phys. Rev. Lett. 78, 2054 (1997), astro-ph/9609169. 
120] M. Zaldarriaga and U. Seljak, Phys. Rev. D55, 1830 (1997), astro-ph/9609170. 

121] M. Kamionkowski, A. Kosowsky, and A. Stebbins, Phys. Rev. Lett. 78, 2058 (1997), astro-ph/9609132. 
122] M. Kamionkowski, A. Kosowsky, and A. Stebbins, Phys. Rev. D55, 7368 (1997), astro-ph/9611125. 
123] W. Hu and M. J. White, Phys. Rev. D56, 596 (1997), astro-ph/9702170. 

124] W. Hu, U. Seljak, M. J. White, and M. Zaldarriaga, Phys. Rev. D57, 3290 (1998), astro-ph/9709066. 
125] U. Seljak and M. Zaldarriaga, Astrophys. J. 469, 437 (1996), astro-ph/9603033. 
126] J. J. Sakurai, Advanced Quantum Mechanics (Addison Wesley, 1967). 
127] R. K. Sachs and A. M. Wolfe, Astrophys. J. 147, 73 (1967). 

128] G. N. Watson, A Treatise on the Bessel Functions. (Cambridge University Press, 1944). 
129] M. Liguori, S. Matarrese, and L. Moscardini, Astrophys. J. 597, 57 (2003), astro-ph/0306248. 
130] E. Komatsu, D. N. Spergel, and B. D. Wandelt, Astrophys. J. 634, 14 (2005), astro-ph/0305189. 
131] M. Liguori et al, Phys. Rev. D76, 105016 (2007), 0708.3786. 
132] F. Eisner and B. D. Wandeh, Astrophys. J. Suppl. 184, 264 (2009), 0909.0009. 

133] N. J. Cornish, D. N. Spergel, G. D. Starkman, and E. Komatsu, Phys. Rev. Lett. 92, 201302 (2004), astro-ph/0310233. 
134] S. Bashinsky and E. Bertschinger, Phys. Rev. Lett. 87, 081301 (2001), astro-ph/0012153. 
135] S. Bashinsky and E. Bertschinger, Phys. Rev. D65, 123008 (2002), astro-ph/0202215. 
136] A. Kosowsky, Annals of Physics 246, 49 (1996), arXiv:astro-ph/9501045. 
137] N. Straumann, Annalen Phys. 15, 701 (2006), hep-ph/0505249. 

138] L. R. Abramo and H. S. Xavier, Phys. Rev. D75, 101302 (2007), astro-ph/0612193. 

139] A. J. S. Hamihon and M. Culhane, Mon. Not. Roy. Astron. Soc. 278, 73 (1996), arXiv:astro-ph/9507021. 
140] A. F. Heavens and A. N. Taylor, Mon. Not. Roy Astron. Soc. 275, 483 (1995), astro-ph/9409027. 
141] A. J. S. Hamilton, (1997), astro-ph/9708102. 
142] E. F. Bunn, Phys. Rev. D73, 123517 (2006), astro-ph/0603271. 
143] M. Kamionkowski and A. Loeb, Phys. Rev. D56, 4511 (1997), astro-ph/9703118. 
144] J. D. Talnian, Special Functions: A Group Theoretic Approach (W. A. Benjamin, 1968). 
145] J. M. Maldacena, JHEP 05, 013 (2003), astro-ph/0210603. 

Appendix A: Integrals of Bessel Functions 
1. Integral of products of spherical Bessel functions 

We will now compute the integrals of spherical Bessel functions that were presented in Section |III| to obtain the 
spacetime window functions. Parts of the methods used here can be found in |28| . The results below also provide the 
motivation for the diagrammatic representation of the window functions shown in Fig. [3j 

The spherical Bessel functions are associated with the matrix elements of the Euclidean group in three dimensions, 
E3 |H|. The Euclidean group E^ consists of the set of transformations that leaves spatial distances invariant - i.e., 
spatial translations and rotations. The rules of group multiplication lead to addition theorems for the special functions 
which realize the group representation, one example of which is the orthogonality condition of Eq. (18 I. Although 
E3 is not compact, the spherical Bessel functions also obey an addition rule, namely |44] : 

(fcr)- ~ ^^ [{kri)ikr2)sm6y- ^n™J(cosyj, ^Aij 

where m is even. 
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2. Integral of three spherical Bessel functions 



We can use Eq. (All and the orthogonality of Legendre polynomials, 

-1 



dxP^"'Hx)Pt'''\x) 



2i+l 



(A2) 



to reduce the product of two Bessel functions to only one Bessel function. Choosing m — 2 due to the demands of 
our particular problem, we have: 



je{kri)je{kr2) = ^ jtikr^) je,{kr2)Se^e> 



t=2 

^ — ^ — / d(cos a) jt {kri)jfj {kr2)Pf'^ (cos a) p/r^' (cos a) 

1 

^ {2e + I) 3 Akn) J I' {kr 2) P^? {cos a) 



d{cosa) P^ '(cos a) (fcri)^ (^'^2)^ sin^ a 



£'=2 



[{kri){kr2) sine 



(A3) 



Now we perform a change of variable, calling r the side of the triangle whose other two sides are ri and r2 , such that 
a is the angle between ri and r2, i. e. r'^ = rf + r2 — 2rir2 cos a. Performing this change of variables, we obtain: 



d(cosa) 



ri+r2 



dr- 



\ri-r2\ rir2 



and therefore: 



ji{kri) je{kr2) = - j drk^ -^-^ J2{kr) P^ ^^(cosa) sin^ a 



Consider, then, the integral that is relevant to us: 



dk ji{kri) jg{kr2) J2{kr3,) = 



rir2 



dr fr 



Pf ' (cos a) sin^ a / dkk"^ J2 {kr) J2 (kr^) 



(A4) 



(A5) 



But now we can employ the orthogonality of Bessel functions, Eq. ( 18 1, so that the k integral gives (7r/2)r S{r — r^) 
and the radial integral can be computed to arrive at the final expresion: 



4 iri,r2,r3) = / dk je{kri) ji{kr2) J2{kr3) = -^^PJ; ^\cosa) sin^ a 

4 Tg 



(A6) 



where ri, r2 and r^ must form a triangle: if they do not, the radial integral yields zero because then r^ cannot be 
equal to some r which, by assumption, forms a triangle together with j'l and r2. Identifying ri — > x, r2 — > Arj and 
rs — >■ Ar/ we obtain the result shown in Eq. (21 1. This integral is also computed in ^28_ . in a more general case but 
employing other methods. 



For the series representation of the window function, consider Eq. (44 1. The same trick that was shown above, i.e.. 



to exchange two spherical Bessel functions for an integral over a Legendre polynomial, can be used to obtain: 
■ jiiki£ri) ji{kar2) J2{kiir3) rir2 f dr 1^2) , ^ . 2 ^ sr^ J2{hir3)J2{kier) 



E- 



km+^{htR) 



dr (_2) , \ 2 1 v^ 

— P^ (cosa)sm a x - >^ 



JJ+iihpR) 



(A7) 



Using now the conjectured orthogonality relation, Eq. (45 I, with I' = 2, we obtain that the window function of Eq. 
(44 1 is indeed identical to the window function of Eq. (21 1. 
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3. Integration of four spherical Bessel functions 

Consider now the integral that appears in the case of A^ = 2 scatterings: 

/•CO 

/f)=/ dkk-^ji{kn)3e{kr2)j2{kr3)j2{kn). (A8) 

Jo 

To benefit from the results obtained above for the case of the integral of three Bessel functions, let ai2 be the angle 



formed by ri and r2, i.e., r = rf + r2 — 2rir'2 cosai2, and use Eq. (Al) to rewrite Eq. (A8l as: 
r(4) 



r ' - 



/•OO 

P^^{cosai2) sin^ai2 dk J2{kr3) J2ikr4) J2{kr) (A9) 

Jo 



rira n^''^ dr 

2 J\ri-r2\ r - ' ' Jo 

fir2 f^^^^^ dr _y o (3) 

' — P^ (COSQ12) sm ai2 X i:^ '(r3,r4,r) 

|ri-r2| ^ 



r(3) 



The integral I2 of three spherical Bessel functions of order two is a symmetric case of Eq. ( A6 ) , and it vanishes 



unless r, r3 and r^ are the sides of a triangle. In this case we can choose any angle in the triangle, so to be consistent 
let's choose that to be the angle between r^ and r4, that is, r^ — fs + ^4 — 2r'3r4COsa34. Since ri, r2 and r must 
also form a triangle, r and the angle 0:34 will be determined in terms of ri, r2, r^, r^ and ai2. The final answer is, 
therefore, given by: 

r w 

4 '(ri,r2,r3,r4) = -rir2r3r4 / ^ P]"^^ (cos 012) sin^ ai2 P2 (cos 0:34) sin^ Q34 (AlO) 

° Jg ^ 

where G is the range of values allowed for r under the conditions that both the triangle of sides (ri,r2,r) and the 
one with sides (r3,r4,r) exist. These conditions are simply the set of inequalities that guarantee that the polygon of 
sides (ri, r2, r3, r4) can exist, i.e., ri < r2 + r3 + 7'4 and the three other cyclical permutations of that inequality. The 
two triangles that must be formed so that Eq. ( A8 1 does not vanish are shown in Fig. p] This result was used in our 



Eq. (281 



4. The general case 

Our problem deals with the propagation of signals between points (events) in spacetime. The first signal propagates 
freely from the source to the point where the photon first scatters, and then there is a set of propagations from one 
scattering to the next, until finally there is a propagation term from the point where the photon have last scattered to 
the point where it is observed. The propagation from the source to the first scattering correspond to a term J2{kArii) 
and each propagation between scatterings to J2{kArii)/{kAr]i)^. The propagation from the last scattering to the 
observation point corresponds to a term ji(kArijq) j {kArjiq)"^ . Besides these propagation terms, the Hankel transforms 
(which were used to go back from Fourier to position space) introduce a ji{kx) into our integral. 

Therefore, for N scatterings we will get integrals over k with an integrand having the following features: 

• two Bessel functions of order £, of arguments kri and kr2 (by convention); 

• N spherical Bessel functions of order 2, of arguments kr^, . . . , krM+2] 

• a factor of /t-2fc-2(JV-2) ^ ^-2(7V-i) 

The method that was used above to compute the integrals in the cases A^ = 1 and N = 2 takes advantage of the fact 
that we can exchange pairs of spherical Bessel functions for Legendre polynomials and radial (or angular) integrals. 
Now, we can do this for every pair of Bessel functions in the TV-scattering integral: if that number is even, every 
Bessel function can be exchanged for an integral over a Legendre polynomial; if that number is odd, an extra Bessel 
function will appear. The final integral over k can then be computed with the help of the lower order integrals. 

For A^ = 3 this procedure leads to: 

/f) = f dkk-^ je{kn) ji{kr2) J2ikr3) j2{kr4) J2ikr5) (All) 
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f dr 12 C dr^i 1^2), x . 2 d(-2)/ n ■ 2 r(3)/ ^ 

/ / Pp ncosai2) sm 012^2 (cos 0^34) sm a^iiyv 12, 'r:^,^^) , 

J ri2 J r34 
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where r^^ makes a triangle together with r^ and r4, and the angles are clearly indicated with respect to their respective 
sides. Notice that, as opposed to the case N = 2, when the sides and the angle ai2 uniquely determines the remaining 
angle of that four-side polygon, in the case A^ = 3 the triangle of sides (r3,r4,r34) is totally free to acquire many 
shapes - see also Fig. [31 It is only when both ai2 and a^^ are given that the angle between ri2 and r^^ is fixed. 
It is also instructive to look at the case N = 4, for the integral of six Bessel functions. In that case we have: 

/f) - Jdkk'^ji{kn)Mkr2)J2{kr3)J2{kr4)J2{kr5)j2ikre) (A12) 

rir2r3r4r5r6 f dr 12 f dr^^^ f dr^Q 



2^ J ri2 J r34 J rse 

X PJ"^'(cosai2) sin^ ai2 P2 {cos au) sin^a34P2 (cosase) sin^ 056 /2 (?'i2,f34,''56) ■ 

The expressions for /j and /) can be obtained in terms of /) ; and so on. With these methods it is trivial to 
compute the spacetime window functions for an arbitrary number of scatterings. 

The set of conditions under which the integrals above are different from zero are those that ensure that each internal 
triangle (corresponding to each instance where two Bessel functions were exchanged for a Legendre polynomial and a 
Bessel function) exist. So, for £ — 4 we would impose: 



ri <r2 + ri2 , r2 < ri2 + ri , ri2 < ri + r2 , 
rs.<ri + r34 , n, < r^i + r^ , r34 < r3 + r4 . 



With the additional condition that ri2 = ^34 in the case £ = A (see the discussion in A. 3), it is trivial to verify that 



these conditions reduce to the inequalities ( 29 1 . 

The set of conditions above simply tells us that the four-sided polygon of Fig. (|3| exists - in other words, that one 
can form a closed polygon with those sides. For any £ the resulting set of conditions ensure that a flat polygon with 
the sides given by ri, . . . , tat exists, i.e.: 

ri <r2+ ... +rN , 

and all cyclical permutations. These inequalities constitute a simple set of constraints that, if not satisfied, imply that 
the spacetime window functions W[^ vanish identically. This simplifies tremendously the integration of the sources 
over time in position space. In particular, one of these inequalities imply that: 

X < Aril + • . • + Ar^jv-i = V - Vn , 

which means that all the sources that contribute to the observables are located at radii x which are inside the past 
light-cone of the observation point at time ry, all the way to the time rj^ when those sources were evaluated, N 
scatterings prior to the observation. 

5. Spacetime window function W^ 

The spacetime window functions regulate how sources at some position x contribute to the observables at time rj. 
For the case of one scattering (order 7), the window function is a (5-function on the PLC, d{x— Ar]), where Ar] = 77 — 77' 
and 77' is the time of the scattering. 

For two scatterings (order 7^), the window function is non-vanishing inside the PLC. In Figs. (p])-(|7| we show a 
few examples of the window functions W^{x, A77; 1). 
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Figure 5: Contour plots of the window functions W({x, Arj; 1) for the cases £ =2, 3, 4 (top panels, left to right panels), and 5, 6 
and 7 (bottom panels.) In these plots x (the radial position of the sources) corresponds to the horizontal axes, and Ar/ = rj — rj' 
corresponds to the vertical axes. For visualization purposes we have fixed A-q' = r] — r\' = \. Physically, this corresponds 
to taking sources at positions x and times 77", and photons which scatter at times 77' before they are observed at time r\. For 
visualization purposes we plotted log |W^/(a;, Ar;; 1)|, so large absolute values of the window functions are indicated by darker 
hues, and the window functions vanish in the white areas. Each lobe corresponds to intercalating negative and positive values 
of the window function. 
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Figure 6: Same as above, but now we fix A77 = 77 — 77' = 1 so the contour plots correspond to Wfix, 1; At;') for the cases 
(. =2, 3, 4 (top panels, left to right panels), and 5, 6 and 7 (bottom panels.) Now A77' corresponds to the vertical axes, and x 
corresponds to the horizontal axes. 
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Figure 7: Same as above, but now we fix a; = 1, so the contour plots correspond to We{l, A77; A77') for the cases £ —2, 3, 4 (top 
panels, left to right panels), and 5, 6 and 7 (bottom panels.) A-q corresponds to the horizontal axes, and Arj' to the vertical 



